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Quasi-classical trajectories study of Ne2Br2(-B) vibrational predissociation: 
Kinetics and product distributions 
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The vibrational predissociation of the Ne2Br2(i?) van der Waals complex has been investigated using the 
quasi-classical trajectory method (QCT), in the range of vibrational levels v' = 16-23. Extensive comparison 
is made with the most recent experimental observations [Pio et at, J. Chem. Phys. 133, 014305 (2010)], 
molecular dynamics with quantum transitions (MDQT) simulations [Miguel et ai, Faraday Discuss. 118, 257 
(2001)], and preliminary results from 24-dimensional Cartesian coupled coherent state (CCCS) calculations. 
A sequential mechanism is found to accurately describe the theoretical dynamical evolution of intermediate 
and final product populations, and both QCT and CCCS provide very good estimates for the dissociation 
lifetimes. The capabilities of QCT in the description of the fragmentation kinetics is analyzed in detail by 
using reduced-dimensionality models of the complexes and concepts from phase-space transport theory. The 
problem of fast decoupling of the different coherent states in CCCS simulations, resulting from the high 
dimensionality of phase space, is tackled using a re-expansion scheme. QCT ro-vibrational product state 
distributions are reported. Due to the weakness of the vdW couplings and the low density of vibrational 
states, QCT predicts a larger than observed propensity for Av' = — 1 and —2 channels for the respective 
dissociation of the first and second Ne atoms. 
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I. INTRODUCTION 

Studying the influence of size-selected solvents on the 
structure and dynamics of molecular systems is essential 
in understanding both molecular energy transfers and the 
transition from the gas to condensed phase. 

In this regard, clusters of rare gas (Rg) atoms doped 
with a diatomic halogen (BC) are particularly convenient 
for at least two reasons: (1) the weakness of the van der 
Waals (vdW) interactions provides a means to effectively 
'separate' the diatom from the environment, which in 
turn allows a relatively simple identification of the dif- 
ferent energy transfer mechanisms at the state-to-state 
level; and (2) the number n of Rg atoms can be spec- 
troscopically selected for their addition induces a known 
blueshift in the vibronic transition B{v') <— X{v = 0)i"— . 

These weakly bound complexes have been the sub- 
ject of intense scrutiny since the pioneering experiments 
of Levy and co-workerai*^!^"— and the theoretical work 
of Beswick, Jortner and Delgado-Barrio^i^ in the late 
seventies. In their experiments. Levy and co-workers 
used laser-induced fluorescence to study Rg„l2 molecules 
based on He, Ne and Ar with up n = 7 I atoms. 

Following the laser-induced vibrational excitation of 
BC, the energy is usually redistributed within the 
molecule leading to the breaking of the vdW bonds. This 
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process is known as vibrational predissociation (VP) and 
provides significant information, for instance, on the dy- 
namics of intramolecular vibrational energy redistribu- 
tion (IVR). Rotational, electronic and other predissocia- 
tion processes are possible. In the former, rotational de- 
excitation of the BC molecule directly provides enough 
energy for fragmentation. In electronic predissociation, 
non-adiabatic transitions to repulsive electronic states of 
BC may instead lead to dissociation of the chemically 
bonded molecule. This was first observed in Arnl-^^ and 
recently studied in HeBr2 and NeBr2^. In general, the 
dependence of the transitions on the details of the po- 
tential energy surfaces (PESs) as well as on the available 
energy makes it possible to extract useful information on 
binding energies and electronic couplings. 

The VP of vdW complexes is by far the predominant 
dissociation process and as such has been the most ex- 
tensively studied. Depending on the vibrational state 
excited, vdW systems fragment following one or several 
elementary steps, which usually include direct dissoci- 
ation, different IVR regimes, evaporative cooling (EC), 
etc. For exainple, the VP of Ne„Br2 clusters is described 
in references [l^ and [III Each step may have different 
accessible final states and characteristic kinetics which 
manifest themselves in experiments through distinctively 
structured product state distributions as in the case of 
AtGI^, HeBr2 and NeBr2i^^. 

All of the above work (and considerably more) has 
proved through the years that, despite their apparent 
simplicity, even small vdW aggregates with n = 1 , 2 un- 
dergo a wide variety of processes that are representative 
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of most of the dynamical pathways observed in more com- 
plex, and conventionally bonded, molecules^. Hence, 
vdW systems have become valuable prototypes in the 
analysis of both inter- and intramolecular energy trans- 
fers. 

For small molecules, exact quantum mechanical calcu- 
lations (EQM) can be performed with modern comput- 
ers and algorithms. Being essentially exact, at least to 
the accuracy of the PESs, these calculations have pro- 
vided a rigorous picture of the VP process. More re- 
cently, they have also addressed in detail the problem of 
IVR dynamics and the role played by continuum reso- 
nances in triatomic systems?'^^-- . Quantum calculations 
on larger complexes almost inevitably use various approx- 
imate methods and there are just a few EQM studies in 
the literature. In 2001, Meier and Manthe studied the 
VP of Ne2l2 using the multiconfiguration time-dependent 
Hartree method^'^ and although vibrational branching ra- 
tios were compared with the experiment^, the main aim 
of their work was to provide a benchmark for future de- 
velopments in the methodology. In 2005, Garcia- Vela 
proposed a full-dimensional, fully-coupled wave packet 
method and used it to study the VP of IIe2Cl2^, ob- 
taining good agreement with experimental lifetimes and 
rotational distributions thus providing a test for the ac- 
curacy of the PES employed. 

Yet, time-dependent, fully detailed investigations of 
energy transfers in vdW molecules containing more than 
one or two Rg atoms still pose a considerable challenge for 
both experiments and EQM calculations. In the former, 
researchers should be able to identify and characterize all 
intermediate complexes, as well as address their dynam- 
ical evolution. In the latter they have to cope with the 
increasing number of degrees of freedom (DOF) and large 
basis sets that eventually make the computational cost 
prohibitive. Thus, in the investigation of larger systems 
hybrid quantum-classical, e.g. molecular dynamics with 
quantum transitions (MDQT) simulationsiSilii^^^, and 
quasi-classical trajectory (QCT) methods have been to 
date the most widely used practical alternatives. 

In fact, most QCT applications have successfully re- 
produced the fragmentation kinetics for at least several 
tri- and tetra-atomic complexes. This is the case for 
molecules such as Ne„l2 (n = 1, 2)^^'^, HelCl^ and 
NeBr2^i^ in which many classical predictions were later 
confirmed by EQM"^^. More recently, some of us have ex- 
plored the extent to which the Gaussian weighted (GW) 
trajectory method^^i^ adds to the applicability of QCT 
in the study of RgBr2 (Rg = He, Ne, Ar) molecules'^^. 
We found that GW may enhance the QCT description of 
product state distributions both quantitative and qual- 
itatively, especially if only a few vibrational states are 
populated or calculations are performed very close to 
a channel closing. However, in cases such as the VP 
of NeCl2^ and generally, when IVR takes place in the 
sparse regime and plays a significant role, it seems that 
only quantum mechanical calculations can be compatible 
with the experimental observations. Hence, despite the 



success of many previous implementations of QCT, the 
validity of classical concepts in the context of the VP of 
vdW aggregates is yet to be clarified. At last, comple- 
mentary studies addressing the dynamical (in)stability 
of these systems exist in the literature. The analysis 
of phase-space bottlenecks in the predissociation of Hel2 
was used as a benchmark in the foundation of phase-space 
transport theory (PSTT)^"— . More recently, irregular 
variations of decay rates as well as details of the absorp- 
tion spectra have been studied through the analysis of 
the phase-space structure and its evolution with increas- 
ing excitation energy^^— . 

In general, testing the accuracy and reliability of theo- 
retical methods largely depends on detailed experimental 
data becoming available. In particular, as stated above, 
realistic kinetic mechanisms can only be obtained if the 
dynamics of all intermediate species is recorded in the 
experiment. This has traditionally proved to be quite 
challenging. For instance, in 1992, Gutmann et al.^^ used 
picosecond pump-probe spectroscopy to study the VP of 
Ne„l2 with n = 2-4 but only registered the formation of 
I2. They therefore inferred the evolution of intermediates 
by fitting the product dihalogen formation to a sequen- 
tial first-order kinetic mechanism. In 2010, Pio et al^ 
reported the characterization of the VP of Ne2Br2(i?) 
at an unprecedented level of detail. Using time- and 
frequency-resolved pump-probe spectroscopy they were 
able to record the real time evolution of all complexes in- 
volved, propose a kinetic mechanism, and determine time 
constants and product vibrational state distributions. 

In this paper, we report QCT and preliminary Carte- 
sian coupled coherent states (CCCS)^i^ calculations on 
the VP of Ne2Br2(S, v' = 16-23). We compare both 
theoretical methods, as well as their ability to reproduce 
the experimental observations^ and previous MDQT 
prediction&ii. Our main goal is to build upon the current 
understanding of the capabilities of QCT in the simula- 
tion of VP processes involving larger vdW clusters, and 
to distinguish whenever possible the results which are in- 
trinsic to the methodology from those that are character- 
istic to the systems under consideration. Complementary 
analysis of full-dimensional and simplified models for the 
VP of the NeBr2 triatomic complex have been very useful 
in the interpretation of QCT results and are conveniently 
discussed. 



II. THEORY 

A. Vibrational predissociation of a tetra-atomic van der 
Waals cluster 

Assuming that the photo-excitation does not provide 
enough energy to break the BC bond, the VP of a tetra- 
atomic vdW aggregate is represented in Scheme [T] where 
the i and f superscripts denote intermediate and final 
states respectively. Two main fragmentation paths are 
possible, which we refer to as the sequential (S) and con- 
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RgiRg2BC(X,«,j) ^ [Rg.Rg^BCiB, v',j')Y 

VP,,, [ Rgi(2) + Rg2(i)BC(B, v\r) "-^^ Rgi + Rg2 + BC(B,^;^/) 

^ S (Ci)Rgi+Rg2+BC(B,«f,/) 

[ (C2) RgiRg2(«L,,7l2) + BC(B,«f,/) 



SCHEME 1: Representation of the VP of a generic RgiRg2BC vdW molecule, showing the (S'i(2)) sequential; as well 
as concerted (Ci) without; and (C2) wit/i molecular formation mechanisms. 



TABLE I. Morse parameters for the various pair interactions 
in Ne2Br2. 







Br2 state 


D/cm-^ 




req/A Ref. 


Br- 


-Br 


X 


24,557.674 


1.588 


2.281 16 






B 


3,788.0 


2.045 


2.667 52 


Br- 


-Ne 


X 


45.0 


1.67 


3.7 20 






B 


42.0 


1.67 


3.9 53 


Ne- 


-Ne 


X, B 


29.36 


2.088 


3.091 29 



certed (C) mechanisms, comprising two channels each. 
The former is associated with two well-defined dissocia- 
tion steps and usually dominates as long as direct disso- 
ciation prevails over IVR. It leads to the loss of two vi- 
brational quanta and have relatively simple implications 
in the VP kinetic o^^i^^ . The latter becomes increasingly 
important as IVR dominates, and is very often linked to 
the loss of more than two vibrational quanta and highly- 
structured product state distributions. All these features 
help in distinguishing between the concerted and sequen- 
tial mechanisms for in the vast majority of cases they lead 
to the same final products (compare 5'i(2) and Ci). Fi- 
nally, when the interaction between the Rg atoms is suf- 
ficiently strong and the dynamical evolution allows for 
favorable configurations, there is a non-negligible proba- 
bility of formation of a RgiRg2 molecule. 



B. Potential energy surfaces 

As seen in Scheme [T] both the ground X and excited 
B electronic states are, at least in principle, involved 
in the VP process. Recent studies on the structure of 
He2Br2 and HelCl complexes^i^ and EQM calculations 
on He2Cl2^^ have shown that the global PES for these 
systems is accurately approximated by 

V = Vagi.BC + VRg2,BC + ^RgiRga + ^BC , (1) 

where Vbc is the interaction potential of the isolated BC 
molecule, VRg.^c (« = 1,2) the vdW PES of the ith 
triatomic aggregate and VRg^Rg^ the potential describing 
the Rgi-Rg2 interaction. 

Here, we extrapolate these results and in a first step, 
we express the global PESs for the Ne2Br2 complex in 



the form ([IJ. It is important to note that, even if this 
proves to be a good approximation for Rg = He and Ne 
due to the weakness of the vdW bonds, neglecting 4- 
body contributions in the global PES of clusters contain- 
ing heavier Rg atoms could eventually fail. Secondly, by 
using pairwise additive potentials for the terms VRg-.BC, 
we neglect 3-body contributions to the triatomic vdW 
PES. Although this is known to be a good approxima- 
tion for the B electronic state, which is where the VP 
process takes place, it is only partially adequate for the 
ground state^i^. In particular, there is theoretical^ 
and experimental^ evidence for an additional minimum 
at linear configurations in the PES of the X state, which 
cannot be reproduced by pairwise interactions. The 
topology of the two electronic states is in fact similar but 
the linear minimum for the B state moves to a longer 
distance and becomes much shallower compared to the 
X state, which makes a pairwise additive description ad- 
equate for the former and not the latter. Nevertheless, 
the calculations presented here refer to the fragmentation 
induced by the photo-excitation of the T-shaped isomer 
which can be correctly reproduced in the pairwise addi- 
tive approximation. All pair interactions are analytically 
modeled by Morse functions, the parameters for which 
have been taken from the literature and are summarized 
in Table m 



C. Quasi-classical trajectories 

1. System model 

A particular Jacobi association diagram yields the 
most convenient set of coordinates to describe the uni- 
molecular dissociation of a triatomic RgBC vdW com- 
plex. Two Jacobi vectors are involved, r which conven- 
tionally runs from the heaviest to the lightest of atoms 
B or C, and R, from the diatom's center of mass to- 
wards the Rg atom. Among its main advantages are the 
symmetrical decomposition of the molecular PES, the ex- 
plicit use of a dissociation coordinate and an associated 
diagonal kinetic operator. It is therefore common when 
studying larger vdW aggregates, i.e. Rgi . . .Rg„BC, to 
choose a 'generalized' set of vectors r, Ri {i — 1, n) which 
inherits most of the aforementioned advantages. These 
are not actually Jacobi co-ordinates but are referred to 
as satellite or bond co-ordinates. Their main drawback 
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have been depicted in Fig. [TJ Choosing J to be zero is 
a well-justified constraint when studying the photodis- 
sociation of rotationally cold species^ as produced in 
the experiment through a supersonic free jet expansion^. 
The complete set of variables in classical phase space T 
is finally obtained by adjusting the respective conjugate 
momenta so that the Hamilton function reads^ 



1 



+Kdw(r,i?i,-R2), 



where 



2^ 



BC 



p2 



Pi 



Ri 



P1P2 



me 



VBc{r) 



mc 
(2) 



(3) 



FIG. 1. (Color online) Coordinate set for a prototypical 
RgiRg2BC vdW molecule in its body-fixed reference frame. 



and 



dof 



Kdw(r, Ri,R2) - V{r, H2) - VBcir) 



(4) 



being a non-diagonal kinetic operator. 

The total angular momentum of a tetra-atomic vdW 
system is given by J = Z1-I-Z2+J, where U is the end-over- 
end orbital angular momentum of atom Rgj with respect 
to BC, and j the rotational angular momentum of BC. 
Using satellite coordinates and keeping J at zero, the 
set of coordinates reduces to {r, Ri, R2,0i,92,(f>}- These 

I 



and 



Here, MRgi,BC = "^Rgi + ("^b + mc) ^ and = -|- 
TO^^ are the inverse of the appropriate reduced masses, 
while the angular momenta 



= ll+ll + 2lrl2 = Pi + Pi + 2Pg, Pg, cos (/) - 2 sin < 



1 



1 



sin^ 6*1 



sin^ 02 



— 2 cos < 



cos Oi COS ^2 



sm t/i sm t)2 

The non-diagonal coupling term. Pi - P2, can be expressed as 



- 2 Pi. 



p2 



Pi 



COS 02 

sin 62 



Pf> 



cos til 
sin Oi 



1,2 



Po. P 



(5) 



(6) 



Pr P2 = PriPr2 cos 7 



R1R2 sin6'i sin 02 
cos 9i cos 02 cos (f) + sin 61 sin 62 
R1R2 

cos 6*1 sin 62 cos cj) — sin 61 cos 62 
R^ 



p2 

Pe^Pe, + 
Pn^Pe, - 



sm 61 cos U2 cos ( 



costal sm 62 



R2 



Pr.Po, 



smyi sin0 smt'TSinffl 

Pr, P4> - ^ . J Pr. P<t> 



i?2 sin 92 



Ri sin 6*1 



cos Oi sm „ „ cos 6^2 sm 

-fg^fcf, „ „ , ^ Pe2Prp, 



R1R2 sin 02 



i?ii?2 sin 6*1 



(7) 



with 

cos 7 = sin 6'i sin 62 cos (/) + cos 9i cos 6*2 . (8) 

Finally, in the interpretation of our results, models 
with two and three DOFs (2/3-DOF) are used for the 
intermediate RgBC complex. The 3-DOF model is a 
full-dimensional J — approximation of RgBC. The 
phase-space variables span a subspace of F, i.e. T3 = 



{r,R,9,Pr,PR,P0}, and all relevant formulae can be 
easily obtained from Eqs. ©-(HI) or the literature, 
e.g. Ref. [s^ The 2-DOF is constructed by additionally 
fixing 9 = corresponding to the equilibrium config- 
uration of NeBr2 in the B electronic state. The phase 
state is thus r2 — {r, R, Pr, Pr} and all necessary for- 
mulae can be obtained by simplification of the 3-DOF. 
In particular, several results can be better understood by 
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analyzing the structure of classical phase space for the 2- 
DOF model. To this end, we have used Poincare surfaces 
of section (SOS), which are a powerful visual tool when 
considering systems with two DOF-—. All SOS employed 
here were constructed from trajectory intersections with 



calculated using 



the hypersurface {r 



,R,Pr > 0,Pr), where rcq is 



the equilibrium bond length of BC. 



2. Initial Conditions 

In general, the QCT simulation of photo-induced pro- 
cesses requires the initial conditions to closely match 
those recreated in the experimenli^. The latter are 
however quantum in nature, and are often quite diffi- 
cult to determine and reproduce classically, especially for 
polyatomic systems. Obviously, the workarounds com- 
monly used in systems with strong inter- /intramolecular 
interactions^ are not applicable here and quantum dis- 
tributions need to be calculated from the molecular wave 
function. 

Although the rigorous form of the initial wave func- 
tion is known to depend on the particular shape of the 
laser pulse, cf. Ref. [20, we assume that the pump laser 
acts during an extremely short time and the system un- 
dergoes a perfect vertical transition. The Schrodinger 
equation for the tetra-atomic system in both the ground 
and excited states is solved using the variational method 
proposed in Ref. H^. Then, the X-state wave function is 
used to calculate the probability distributions (statistical 
weights) for all relevant variables, while the eigenvalues 
from the B-state determine the VP energetics, as in our 
previous work^i*^. Once the X- and i?-state quantum 
problems are solved, QCT initial conditions are calcu- 
lated following the algorithm proposed in Ref. [6^ (with 
only a few minor modifications): 



1. a set of 'quantum' numbers {^i, j} is selected us- 
ing their respective distributions — checking it satis- 
fies the triangle condition A(/i, l2,j) to ensure com- 
patibility with J — 0; 

2. the angles 9i, 62, 7 are generated following the 
appropriate distributions, and (j) calculated from 
Eq. ®; 

3. the system of non-linear equations and © is 
solved for the angular momenta Pg-^, Pg^ and P^; 

4. r is generated from the calculated distribution, and 
the value of j used to calculate Pr from Eq. ([3]). 
Here, Hbc is conveniently replaced with , 
i.e. the eigenvalue corresponding to the (i?, v' ,j' — 
0) state of BC^^, and the sign of Pr randomly cho- 
sen; 

5. i?i and R2 are generated with the corresponding 
distribution, and approximate conjugate momenta 



'2^r 



Rgi.BC I "^Rgi^C 



B.v' 



2/i 



Rgi,BC 



— VR,gj,BC 



where £, 



Rgi,BC 



(9) 

is the eigenvalue corresponding to 
the RgiBTJ triatomic aggregate with BC in its 
{B, v') stateS? , and the sign of the momentum is 
once again randomly chosen; 

6. finally, the energy associated to the non-diagonal 
kinetic term P1P2 and Vng^Rgj is redistributed 
within the vdW modes by writing Pr^^^) ~ 
Pr'^'2 ^1(2)' where Ai and A2 are quan- 
tities to be determined. Their calculation 
requires evaluating 7i'^^° = , the en- 

ergy of the tetra-atomic complex with Br2 
in the {B,v') state, using the set of values 
{r, iii, i?2, ^1, 02, 0, Pr, Pflf , Pflf , Pe,,Pe, , Pj.}: as 
well as the additional constrain P^^/P^^ = 
A1/A2. 

Steps 1-6 are repeated until A'tot initial conditions are 
obtained. 



3. From classical magnitudes to observables 

We have recently discussed the statistical handling of 
QCT results for comparison with experimental observ- 
ables in the specific case of the VP of triatomic vdW 
systemSi^. There is however one fundamental difference 
between the process in tri- and tetra-atomic molecules, 
namely, the possible formation in the latter of an interme- 
diate complex. For consistency, the general methodology 
is only summarized here with emphasis on the changes 
made in order to analyze the QCT results in our specific 
case. 

When the molecule completely dissociates, Kdw ^ 
and both j and Hbc become integrals of motion with the 
latter corresponding to the classical energy of the BC 
fragment, Ebc- The final semi-classical ro- vibrational 
quantum numbers for the ith dissociated trajectory, 
i.e. j and j^j, are then rigorously constant and are 
given by 

f 1 / , 1 
2TTh T 2 



V2^BC 



BC, 



V,s dr - -, (10) 



where Vcs is the effective interaction potential of BC (in- 
cluding the centrifugal term) and the closed integral is 
evaluated over one BC vibrational period^ (''min and 
Tmax are the classical turning points), while 



Jc.i 




(11) 
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which results from jf = fc_i{jl i + l)h^^^. Due to the 
weakness of the vdW interaction, the BC vibrational 
state in the intermediate complex, wj, j, can be estimated 
using Eq. (ITUl) . This requires replacing the final state 
magnitudes with approximate values for the intermedi- 
ate state, despite these not being rigorously defined. 

In what follows, we neglect the quantized nature of the 
rotational DOF as it has a very high density of states and 
refer simply to channel Av' = — v' by the associated 
final vibrational quantum number, . 

In the procedure known as histogram or standard bin- 
ning (SB), the probability density P{M,v^) that a given 
observable has the value M and the final BC vibrational 
state is can be formally written as 

PsB(M,«f) = J dTp{T)d[M{T)-M]E[v\Ty,v',l], 

, (12) 

where p(T) is the probability distribution of the initial 
phase space state, F, and 

E{x;x*,A) = (a;* 4- A/2 - x) e (x - a;* + A/2) . 

(13) 

S{x) and Q{x) are respectively the Dirac and Heaviside 
functions. M{T) and v^{T) are the final values of the 
observable and the vibrational action in terms of F. It is 
relatively easy to see from Eq. ([T5)) that S defines a square 
barrier function of x, which equals 1/Aon [x* ~A/2,x* + 
A/2] and everywhere else. Also, Psb in Eq. is 
normalized to 1. 

In practice, Psb is estimated by means of the Monte- 
Carlo expression 
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S[Af(F,);Mfe,afc] 

Z— 1 

xE[v'{T,);v',l], (14) 



where Mk is the kth midpoint in the {Afj^j^ set par- 
tition of the interval [Mmin, Mmax], afe = Mk+i - Mk 
and A'diss is the total number of dissociated trajectories, 
which is assumed to be large. 

In the GW procedure one simply replaces the square 
barrier function S [ti^(F); , l] in Eq. with the Gaus- 
sian 

9[vHr);v^,e]^ ^ L, (15) 

in which e is usually kept at 0.05. Hence, Pqw is simply 

PGw{Mk,v') « -— V S[M(F,);Mfe,afc] 



<g[/(F,);^;^e] 



(16) 



In the tetra-atomic case, variations of these general ex- 
pressions can describe either the intermediate or final 



states. To do this, it is sufficient to use the appropri- 
ate weight function for the desired state for the sec- 
ond S function in Eq. ([T4l) and the Gaussian function 
in Eq. p3|) respectively. We must note that the weight 
for the intermediate state is calculated as the product 
of the weights corresponding to the semi-classical vibra- 
tional level of BC and that of the intermediate vdW com- 
plex: w' = w\ (uj.) W2 , where w refers to either S or 
g. This could dramatically increase the number of tra- 
jectories necessary for convergence. Moreover, although 
Eq. (fTU)) provides a simple means to estimate the semi- 
classical analogue for the vibrational state of BC, that 
of the vdW complex is much more difficult to evaluate. 
A simpler approach, known as the 1GB procedure, uses 
the total ro-vibrational energy rather than that of the 
individual states to calculate the corresponding weight: 

w' = w' ^£^^g"Bc) ■ Both the and S^g^^c eigenval- 

ues are involved in this calculation, see Sec. Ill C 21 The 
method was proposed by Czako and Bowman^^ on the 
basis of rather intuitive arguments and was later theo- 
retically validated by Bonnet and Espinosa-Garcfa^. In 
the particular case when the first and second dissocia- 
tion steps are statistically independent, the total weight 
can be calculated as the product w = w^w^ . In addition, 
mixed strategies may be tested by simply using w = w'*^^^ , 
hence calculating the desired state distributions indepen- 
dently of the final/intermediate vibrational state. 

The time evolution of the populations of the various 
complexes can be obtained by small modifications to the 
formulae determining survival probabilities, i.e. the prob- 
ability that a given complex has not dissociated at time 
t, which read^ 



PsB{t) - 1 



1 



JVdi. 



P( 



GW 



(0 = l-=^lv^E3.e(i-tO, (17) 
2^1=1 9i i=i 



with gi = X^vt) ^' ^) ^"^^ dissociation time 

for the zth trajectory. In the particular case of a Rg2BC 
cluster, for which ti and t2 are the dissociation times for 
the first and second Rg atoms, the different populations 
are exphcitly obtained as: (1) P4{t) — P{ti), the decay 
of the parent molecule; (2) P3{t) — P{t2) — P{ti): time 
evolution of the intermediate RgBC population; and (3) 
P2{t) = 1 — P{t2), formation of the BC fragment. More- 
over, (4) P^{t) = P{t2 — ti) directly represents the decay 
of the RgBC complex. The reasons for the subscripts will 
become apparent in Section IlII Al 



4. Simulation details 

Batches of A'tot = 5 x 10^ trajectories were propagated 
using an adaptive-stepsize Bulirsch-Stoer method^S. 
Each trajectory was followed until one of two condi- 
tions was fulfilled: (1) the two Ne atoms dissociated. 
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i.e. Ri{2) > ^diss"^ — 14 A; or (2) the propagation time 
t = Tniax = 1200 ps. Dissociation is said to occur when 
Ri(2) — 14, the distance at which the vdW interactions 
become neghgible. The integration parameters were ad- 
justed so that the maximum error in total energy did not 
exceed AE = 10^^ cm"\ i.e. less than 10^^ E. With 
this choice of parameters all trajectories dissociated at 
least one Ne atom while over 88-99%, depending on v', 
completely fragmented. A typical trajectory requires less 
than 1 second of CPU time on an Intel® Core^" 17 Q720 
(6M Cache, 1.60 GHz) processor. 



D. Cartesian coupled coherent states 

The CCCS method is a trajectory-based quantum dy- 
namics technique designed to be similar to classical tra- 
jectory simulations. The main focus of this paper is the 
QCT predictions and thus CCCS results are mainly pre- 
sented to corroborate the conclusions and the quality of 
the QCT results. Consequently, we will only give a brief 
outline of CCCS in order to aid comprehension of the re- 
sults that are included. The interested reader is referred 
to Ref . which not only describes in detail our previous 
work on the VP of NeBr2 but also outlines the exten- 
sion to larger clusters. A more detailed CCCS study of 
Ne2Br2 will be published in due course^. 

The QCT method described hitherto is concerned with 
the point-like nuclei that make up the Ne2Br2 cluster 
whereas the CCCS method allows us to study the time 
evolution of the associated quantum mechanical wave 
function. The CCCS method expands the wave function 
using a basis of coherent states (CS). In CCCS, the CS 
are Gaussian-shaped wave packets that describe both the 
position and momentum for each Cartesian DOF of each 
atom. Thus for Ne2Br2, each CS has 24 dimensions. The 
CS move on the same PES, and according to the same 
equations as the nuclei in QCT. However, the potential 
energy for the CS is the convolution of the CS and the 
PES. In contrast, within the QCT method, the potential 
energy of a given configuration is simply the value of the 
PES at that point in configuration space. In CCCS, the 
convolution may be done in advance of the simulations 
and it gives the so-called averaged or re-ordered PES (or 
equivalently, Hamiltonian) upon which the centers of the 
CS move. 

The amplitudes evolve with time according to an ex- 
pression derived from Schrodinger's equation and which 
depends upon the fact that the basis functions overlap. 
This overlap couples the amplitudes of the basis functions 
together. If the basis functions are separated by a suffi- 
ciently large distance in the 24-dimensional phase space, 
they cease to overlap and become decoupled. Once all 
the basis functions have decoupled, the CCCS method 
is essentially a semi-classical method. However, even in 
its semi-classical limit, the CCCS technique takes into 
account the majority of zero point energy effects. 



1. Simulation details 



The initial coordinates for the basis functions are de- 
termined from a set of bond lengths and momenta so that 
the total linear and angular momenta of the cluster are 
zero. For Br-Br, these were chosen from the phase-space 
trajectory at the energy expectation value of the isolated 
molecule. In the case of Br-Ne and Ne-Ne, they were 
chosen randomly from the ground state wave functions 
of the bonds. Initially, each basis consisted of 400 ba- 
sis functions and the results have been averaged over 10 
different basis sets per Br2 vibrational level. 

As CCCS simulations are considerably more expensive 
than QCT calculations, a Ne atom is said to have dis- 
sociated once the mean distance between it and the Br2 
molecule exceeds -R^isfP^ ~ 10 A. Similarly, the second 
stage dissociation is said to occur when the distance be- 
tween Br2 and the remaining Ne exceeds the same cutoff. 

In our previous work^— , we found that the basis func- 
tions in the simulations of the VP of NeBr2 remained 
coupled until the clusters completely dissociate. In the 
current work, however, the basis functions decouple much 
more quickly due to the increased dimension of phase 
space. Obviously, the time at which each basis function 
decouples varies greatly but it is typically less than 10 ps: 
the mean time in one simulation was about 6 ps. As will 
be discussed later, this leaves comparatively few inde- 
pendent basis functions to describe the dissociating wave 
function and leads to increased lifetimes. 

The calculated lifetimes can be improved by spawning 
a secondary simulation each time a basis function perma- 
nently decouples in which the decoupled basis function 
is expanded on a new basis of 50 CS. The calculation 
of the lifetimes is thus a two-part process: first we per- 
form a simulation of the whole wave function expanded 
on the basis of 400 CS; and then we perform 400 sec- 
ondary simulations each with 50 basis functions describ- 
ing one of the basis functions of the original simulation 
starting from when that basis function decoupled. We 
found that decoupling invariably occurs before dissoci- 
ation which makes the combination of the dissociation 
curves from 400 secondary simulations easier. 

The time required to perform the CCCS calculations 
depends upon the length of time for which the basis func- 
tions remain coupled and increases nonlinearly with the 
number of basis functions. As an example, the mean 
time for v = 17 for 100 ps (10^ steps), using 400 ba- 
sis functions was 3.2 days. The mean time for the 
re-expansions however was 1.3 hours per simulation on 
Intel®Xeon® 'Woodcrest' processors. 
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III. RESULTS AND DISCUSSION 



written as 



A. Kinetics, lifetimes 
1. QCT 

PSTT establishes that certain structures in classi- 
cal phase space (resonant islands, tori, cantori. . . ) act 
as inter- and intramolecular bottlenecks to the diffu- 
sion of trajectories from regions defining 'reactants' and 
'products— A rigorous (classical) kinetic mechanism 
must therefore reflect the details of these partial ob- 
stacles. In addition, a complete rationalization of the 
energy-dependence of decay rates requires the stability 
of periodic orbits and the evolution of bifurcations in the 
phase-space portrait to be addressed in depth. While 
this is possible and has already been carried out for sim- 
ple 2-DOF models of molecular systems, including vdW 
complexes^i^i^Sr— , its extension to more DOF intro- 
duces serious technical difficulties^ and constitutes in 
itself an active field of research. In fact, further below we 
use some of these results and a 2-DOF model of NeBr2 
to better understand the relation between product states 
and dissociation lifetimes, as well as certain features in 
the time-evolution of the intermediate complex. 

Analyzing the phase-space structure of Ne2Br2, or even 
NeBr2, is however far from our objective here and, at this 
point, we take advantage of one important experimental 
conclusion: At least for low v' levels, a simple sequential 
kinetic mechanism (Section III A[) provides an excellent 
fit to the delay scans (Figs. 7 and 8 in Ref. j47)). This 
is consistent with several results from our QCT calcula- 
tions: (1) once the Ne2Br2 complex decayed into NeBr2, 
the only process observed was the fragmentation of the 
resulting triatomic into Ne and Br2 products; (2) less 
than about 0.05% of all trajectories dissociated via the 
concerted mechanism (subset of all trajectories for which 
both Ne atoms dissociated within one Br2 vibrational pe- 
riod, excluding those predicted by the sequential mecha- 
nism to dissociate within that time interval^^); and (3) no 
significant statistical correlation was found between the 
variables defining the intermediate complex and ^2 — ii- 
The latter confirms, at least from the kinetics viewpoint, 
that the second dissociation step is nearly statistically in- 
dependent on the first^i. All of these observations agree 
with the fact that: (4) EQM calculations on the VP of 
NeBr2 show that direct predissociation dominates within 
the vibrational range explored here^S. 

The proposed kinetic mechanism is schematically 



7V4 



No 



(18) 



where N4, N3 and N2 stand for the Ne2Br2, NeBr2 and 
Br2 populations respectively, while ri and T2 are the cor- 
responding lifetimes. Mathematically, Eq. ((T5)) can be 



with ki ~ Tj 

N2{t) 



dN^/dt = -kiNi, 
diVa/dt = fciiV4 ~ k2m, 

dN2/dt = k2N3, 

and its solution 



(19) 



iV4(0)exp(-fcii), 
A:iiV4(0) 



Ntot- 



k2 - ki 



exp (— fc2<) + 



k2 - ki'' 



(20) 



where iV4(0) = iVtot and 7V3(0) = N2{0) = 0. Given the 
dependence on ri and T2 of Eqs. (f20| . various schemes 
to extract the lifetimes from survival probability curves 
become possible. One can: (1, 2) first obtain ti from 
fits to Ne2Br2 decay curves and use this value to get 
T2 from the time dependence of NeBr2, or alternatively, 
Br2 populations; or (3, 4) directly infer both values from 
the evolution of NeBr2 or Br2 populations. Alternatively, 
taking into account that the individual dissociation times 
ti and t2 are nearly uncorrelated, one may: (5) directly 
extract the value of T2 by fitting the probability of com- 
plete dissociation after a time t2i = t2—ti (fragmentation 
of the triatomic once it is formed) to a single exponen- 
tial decay. This is implicit in our kinetic mechanism and 
readily seen by adding the first two equations in (|19l) . 
Schemes (l)-(4) are somewhat similar to the experimen- 
tal freedom of determining the lifetimes by measuring the 
time dependence of the Ne2Br2/NeBr2 signal in a given 
Br2 vibrational state, referred to as 'disappearance' of 
the parent complex, and that of the NeBr2/Br2 signal 
in the vibrational states resulting from dissociation, the 
'appearance' of the daughter molecule, etc. Depending 
on the signal monitored, different values are obtained in 
the experiments (and EQM), which provide a measure 
of IVR^'^'^s. However, QCT will yield the same life- 
times no matter what scheme is used (though numerical 
fitting and convergence will in practice result in small dif- 
ferences). This is easily understood for according to our 
criteria, the disappearance of Ne2Br2 exactly matches the 
appearance of NeBr2, etc. As an example, the relevant 
probability curves, i.e. Pi{t) = Ni{t)/Ntot, i = 4-2 and 
P^{t), corresponding to v' = 21 are shown in panel (a) 
of Fig. [21 together with the respective fits. We have ver- 
ified that the various schemes yield similar results, and 
experimental and theoretical lifetimes are summarized in 
Table |TT1 In the latter, column 'w" corresponds to the 
vibrational level in the B electronic state to which the 
Ne2Br2 molecule is excited by the laser pulse. Column 
'Species(D)' then corresponds to the molecular product 
(in the specific vibrational state v) that is monitored in 
the experiment. 'ti' and 'r2' are the lifetimes of Ne2Br2 
and NeBr2. In particular, as seen in Eqs. (I^H)) . the frag- 
mentation of Ne2Br2 depends only on fci, which is why 
no T2 is reported in the experiment for this molecule. 
As discussed above, the theoretical methods yield nearly 
identical results independently on the probability curve 
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FIG. 2. (Color online) Time evolution of the population of all complexes involved in the vibrational predissociation of Ne2Br2(-B, 
v' = 21): (a) QCT; and (b) CCCS calculations. 



used for the fitting, and a single pair of values n , T2 is 
thus reported for the system (note that no T2 value was 
reported in Ref. [TTI ) . 

Simple inspection of panel (a) in Fig. [2] shows that 
the agreement between our fits and QCT calculations is 
very good, which confirms the accuracy of the kinetic 
mechanism proposed. This is the general trend for all 
vibrational levels explored here. Even so, in most cases, 
small discrepancies occur at short and large times. From 



TABLE II. Experimental-i and theoretical — this work and 
MDQTii — lifetimes (in ps) in the vibrational predissociation 
of Ne2Br2(B,w' = 16-23). 



Experiment 


QCT 


CCCS 


MDQT 


v' Species(«) n 


T2 


Tl 


T2 


Tl 


r2 


Tl 


16 Ne2Br2(16) 




62.1 


134.0 








17 Ne2Br2(17) 32±3 




48.5 


115.9 


44.0 


98.9 


29 


NeBr2(16) 30±3 


88±3 












Br2(15) 31±2 


82±3 












18 Ne2Br2(18) 28±3 




33.8 


89.6 








NeBr2(17) 27±5 


58±5 












Br2(16) 28±3 


55±4 












19 Ne2Br2(19) 




28.0 


66.5 






24 


20 Ne2Br2(20) 




23.9 


51.1 






19 


21 Ne2Br2(21) 16±3 




19.3 


40.5 


21.7 


53.1 


16 


NeBr2(20) 14±5 


30±2 












NeBr2(19) 17±4 


54±3 












Br2(19) 16±2 


29±2 












Br2(18) 16±2 


47±2 












22 Ne2Br2(22) 




16.2 


34.5 






13 


23 Ne2Br2(23) 




14.1 


28.9 


15.3 


39.4 





the comparison of the classical survival probability curves 
corresponding to the decay of the Ne2Br2, red curve, 
and the fragmentation of the intermediate NeBr2, green 
curve, the nature of the plateau usually observed at short 
times becomes rather clear. It is simply an artifact re- 
sulting from both the classical description of the process 
and the way initial conditions are sampled. More explic- 
itly: there is a minimum time for the gradual classical 
energy transfer from the Br2 vibrational mode to disso- 
ciate one Ne atom which, in the case of the first decay, 
has been located close to the vdW minima at the belt-like 
configuration. In contrast, a quantum vibrational tran- 
sition may immediately release enough energy for disso- 
ciation to occur, as readily seen in CCCS curves from 
panel (b). The fact that no such feature is observed in 
the second QCT step further confirms this reasoning for 
once the first atom is lost, the dynamics would sample 
more evenly the available phase space and have placed 
the second Ne atom arbitrarily close to i?diss- One may 
say that this apparent non-exponential behavior (in the 
fragmentation kinetics) mainly arises from initial state 
selection^^ . 

The reasons for the discrepancies at larger times are 
more complex and better understood using a 2-DOF 
model for the VP of NeBr2 (with 6* = tt/2). Analysis of 
the SOS shows a 1:10 nonlinear resonance (for v' = 21) 
surrounding the stable central region around the point 
(i? = Rcq, Pr = 0). This resonance occupy a signifi- 
cant proportion of the available phase space. The extent 
to which such intramolecular bottlenecks affect the over- 
all time evolution of the classical ensemble depends on 
the proportion of initial conditions lying inside or rel- 
atively close to their boundaries. Actually, early work 
on PSTT showed that kinetic mechanisms can be con- 
veniently modified to account for such behavior^. In 
our case, only a relatively low percent of trajectories, 
about 0.3% for v' = 21 in Fig. [51 show a strong non- 



10 



exponential behavior, always after the dissociation of 
the first Ne atom. This is most likely due to remain- 
ing lower-dimensionality tori, which have been recently 
demonstrated to play a significant role as bottlenecks be- 
tween diffusive and statistical behavior in systems with 
more than two DOF— . 

Given all of the above, the actual equations used to 
produce the fits shown in Fig. [5] are in fact slightly 
modified versions of Eqs. (1^01) which account for both 
sources of non-exponential behavior. That at low times 
is 'avoided' by fitting from > 0, while the long-time 
behavior is modeled by adding a given constant (which 
equals the asymptotic proportion of non-dissociated tra- 
jectories at each v'). We should note that, after adding 
this constant, all coefficients need to be adequately mod- 
ified to recover the correct behavior at t = 0. 



2. CCCS 

The survival probabilities for the different complexes 
for v' = 21, as calculated using CCCS, are shown in 
panel (b) of Fig. [21 In contrast to panel (a): (1) dotted 
lines (labeled '400 CCS') correspond to the initial CCS 
simulations with 400 basis functions, the re-expansion of 
which gives the main curves; (2) the survival probability 
of NeBr2 once it is formed from the Ne2Br2 cluster is 
not presented, as this time is ill-defined in the CCCS 
calculations; and (3) dot-dashed lines (labeled W^^^ — 
16 ps') uses the experimental valued to model the time 
evolution at short times. The fits are obtained using the 
same sequential kinetic mechanism as for QCT, and the 
lifetimes given in Table |lT] are from fitting the modified 
versions of Eqs. (l20l) to the CCCS data, also necessary 
here, mainly to account for those basis functions that do 
not dissociate. 

Theoretically, CCCS will provide a quantum mechani- 
cal description as long as the basis functions remain cou- 
pled, which ideally should be for a time comparable to the 
process of interest. As in our previous study of the NeBr2 
system^, CCCS removes the nonphysical plateau shown 
by QCT at short times, thus accurately reproducing the 
initial quantum dynamics. However, it becomes harder 
to fit the results of CCCS to a simple model of sequential 
dissociation at longer times. This could be because quan- 
tum mechanics makes non-sequential dissociation more 
likely. Previously, non-sequential IVR-EC contributions 
to dissociation were found to be important in MDQT sim- 
ulations for V > lA where they account for one fifth of all 
dissociative trajectories^i. Also, experiments^^ seem to 
indicate that complicated non-sequential mechanisms are 
important for v > 19. However, it would be premature 
to claim that the current CCCS simulations prove the 
pertinence of non-sequential mechanisms. Although our 
re-expansion technique (Section [II D 1|) greatly increases 
the time scale at which CCCS works, it is still below the 
range of hundreds of picoseconds required for the descrip- 
tion of the whole fragmentation process. At such long 



time scales, CCCS works as a semi-classical technique 
running largely uncoupled CS on the averaged potential. 
It accounts for zero point energy effects but not for the 
full quantum coupling in phase space. We are currently 
working on strategies to further extend the coupling time 
between the CS in the quantum simulations. More de- 
tailed (and costly) investigations are necessary and the 
results will be published elsewhereSS. 

In general, both QCT and CCCS predictions are in 
very good agreement with each other, thus CCCS con- 
firming the good quality of much simpler and less expen- 
sive QCT results. Both methods are also in good agree- 
ment with the experimental lifetimes, over the whole 
range of initial vibrational states. However, although 
providing a near-quantitative description of the kinet- 
ics, they seem to consistently yield larger than observed 
values. In the case of QCT, this may be intrinsic to the 
methodology and the existence of low-dimensional tori in 
classical phase space. Additional CCCS calculations are 
needed in order to verify whether this is inherent to the 
method or due to, e.g. an inadequate representation of 
the PES or the wave function. 



B. PSTT, SB, GW and QCT dissociation times 

Based on PSTT, it seems reasonable to expect that tra- 
jectories with the highest Gaussian weights, which con- 
tribute the most to the GW curves and correspond to 
vibrational actions the closest to integer values, do not 
necessarily explore the same regions in phase space than 
the whole ensemble and could therefore undergo different 
kinetics. Hence, at least in principle, applying the SB or 
GW methods can produce different time evolution prob- 
ability curves (Section III C 31) and correspondingly, result 
in different lifetimes. In turn, GW results should be the 
closest to the experimental values for this method effec- 
tively eliminates trajectories which finish in nonphysical 
vibrational channels. Quite surprisingly, no significant 
differences between SB and GW lifetimes were found in 
our previous calculations on the VP of RgBr2 (Rg = He, 
Ne, Ar) molecules (see Fig. 1 in Ref . [35l) . where several 
possible reasons were proposed. 

In the VP of Ne2Br2 studied here, the predictive capa- 
bilities of GW can be additionally tested for the experi- 
ment provides detailed information on the kinetics associ- 
ated to different vibrational channels for v' = 21 (see Ta- 
ble [Hi . Unfortunately, once again, our calculations seem 
to corroborate that both the SB and GW methods result 
in nearly identical lifetimes, with no significant variation 
between the various vibrational channels. In other words: 
the GW procedure was found to be unable to provide vi- 
brational resolution in the dissociation lifetimes. In the 
rest of this section, we will show that 2- and 3-DOF mod- 
els for the VP of NeBr2 are very helpful in the detailed 
analysis of this issue. 

The relation between trajectory lifetimes and final 
states is rather intricate, as is apparent from Fig.[31 Panel 
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Decreasing dissociation time 



FIG. 3. (Color online) Dependence between dissociation life- 
times and final vibrational states in the vibrational predisso- 
ciation of NeBr2 from v' = 23: (a) 2-DOF; (b) 3-DOF model. 



(a) shows how final vibrational states relate to the indi- 
vidual dissociation times for an ensemble of 5 x 10'* tra- 
jectories in the 2-DOF model. As readily seen, there is 
a relative insensitivity between these magnitudes for tra- 
jectories with longer to medium lifetimes, while a spe- 
cific pattern arises for shorter dissociation times. These 
finger-like structures can be understood by analyzing the 
VP process in phase space, in the spirit of PSTT. For sim- 
plicity we will only introduce here the essential elements 
of the methodology that are required by our interpreta- 
tion. Further information can be found in the seminal 
work of Davis et ai, most notably in Ref. 37. 

Figure S] depicts the main structures in the SOS of a 
2-DOF model for a generic triatomic. The SOS is con- 
structed by the procedure described in Section III CI and 
the molecule is assumed to dissociate if i? > i?diss- The 
separatrix (continuous red line) encloses the interaction 
region and constitutes the intramolecular bottleneck to 
dissociation. The fingers (green/blue dashed lines) limit 
the region available to the outgoing/incoming flux. All 
dissociating trajectories will intersect the region enclosed 
by the outgoing fingers (dashed, green) and the separa- 
trix, once only. Thus, after leaving the interaction region, 
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FIG. 4. (Color online) Schematic view of the main struc- 
tures in the SOS of a 2-DOF model for the VP of a generic 
triatomic: separatrix (red, continuous line), outgoing (green, 
dashed line) and incoming (blue, dashed line) fluxes. 



any given trajectory will inevitably dissociate and con- 
secutive intersections will lie in different fingers as i?(i) 
increases monotonically. All that rests now is to notice 
that trajectories closer to the fingers' 'tips', i.e. those 
with the largest Pr, correspond to the smallest v avail- 
able (and vice versa). The previous follows from the fact 
that the total energy E is fixed in the ensemble. At last, 
the existence of finger-like structures like those in Fig. [3] 
will strongly depend on the characteristics of the dynam- 
ics within the interaction region. If completely stochas- 
tic, initial conditions will be effectively 'forgotten' and 
trajectories will access the outgoing fingers after explor- 
ing the interaction region during a random propagation 
time. No correlation between dissociation times and final 
vibrational states will be observed. On the other hand, if 
intermolecular couplings are not strong enough, a subset 
of trajectories whose initial conditions lie within a spe- 
cific region inside the separatrix will directly accesses the 
outgoing fingers. These trajectories will therefore disso- 
ciate rapidly and a pattern will arise. Following this line 
of reasoning, it is not difficult to see why the structures 
in panel (a) of Fig. [3] have precisely their shape. Panel 
(b) in the same figure shows results for the 3-DOF model 
of NeBr2. The similarities between the 2- and 3-DOF 
models are evident, although the finger-like structures 
get increasingly blurred the larger the number of DOF. 
The number of such finger-like structures depends on the 
kinetic energy available to dissociation, hence decreasing 
with increasing vibrational excitation. Additional cal- 
culations show that all structures completely disappear 
for energies above the linearization threshold (Br-Ne-Br 
configuration barrier). 
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TABLE III. Energies (in cm ^) in the vibrational predissoci- 
ation of Ne2Br2(-B, f ') (see text for full details). 



v' 


A c.S,d'-1 




j-)B,v' 


rpv —1 
avail 


rpv —2 
avail 


16 


107.4 


63.3 


137 


33 


81 


17 


103.8 


63.1 


137 


29 


74 


18 


100.2 


62.8 


137 


26 


67 


19 


96.6 


62.6 


137 


22 


60 


20 


93.1 


62.3 


136 


19 


54 


21 


89.5 


62.0 


136 


15 


47 


22 


85.9 


62.5 


135 


13 


40 


23 


82.3 


62.4 


135 


9 


33 



C. Product state distributions 



Ro-vibrational product state distributions were calcu- 
lated from the ensemble of quasi-classical trajectories as 
described in Section III C 31 As many relevant features 
of these observables strongly depend on the energetics of 
the VP process, details of the latter are summarized in 
Table IIIII (1) the vibrational energy gap associated to 
the Av' = -I channel, A£^;^ = - (2, 

3) NeBr?^ and Ne2Br2 binding energies, jj^'^^''^'^^'^ ) — 

~^NcBr2 ^^'^ ^0 ^ — ; and (4, 5) the energies 

available after dissociation of the first and second Ne 
atoms via, respectively, the Av' — —1, —2 channels, 
^.e. Ei-J and E^. 

Tests on basis-set convergence yield an estimate for 
the accuracy of our D^'^ values of about 1 cm~^ (any 
further improvement was considered unnecessary for our 
purposes here). These theoretical values are just slightly 
different from the 141 cm~^ estimate of Pio et ali^, 
and suggest an interaction energy for the Ne-Ne bond 
of about 10-12 cm~^, which is to be compared with the 
17 cm~ of isolated Ne2 used in Ref. Our predictions 
are, however, consistent with the value of 12.81 cm^^ in 
Ref. [ill for the effective Ne-Ne bond. Moreover, our cal- 
culations predict an average Ne-Ne interatomic distance 
at the minimum energy structure of about 3.2 A (within 
3.8% of the value for Ne2, see TableU]), also in good agree- 
ment with the estimate in Ref. It is important to note 
that all values in Table Hill except for ^^°B'2(s,i> )^ 
theoretical predictions based on the B-state PES used 
in this work. This is the reason for small discrepancies 
with some experimental predictions'*^. For instance, we 
estimate the closing of the Av' = — 1 channel for the 
dissociation of the first Ne atom to occur at v' = 25, 
as opposed to the experimental v' = 23. Also, complete 
dissociation via the Av' = — 2 channel would be energet- 
ically accessible up to v' = 27 (instead of v' = 25). 



1. Vibrational branching ratios 

Experimental and QCT vibrational branching ratios 
after the dissociation of one Ne atom, NeBr2(u' — 
n):NeBr2(w' — 1), and two Ne atoms, Bi2{v' — n):Br2{v' — 
2), are given in Table HVl It is readily apparent that QCT 
fails to capture the physics behind the vibrational distri- 
butions in the VP process, irrespectively of the statisti- 
cal procedure (SB or GW) employed. This is a general 
effect of the large mismatch in the strengths of the dif- 
ferent DOF involved in the VP of vdW clusters, at least 
for moderate vibrational excitations of the chemically- 
bounded diatom. The same behavior may therefore be 
expected in similar systems. In the absence of resonances, 
such a mismatch causes the energy transfer from the vi- 
brational to the vdW dissociation modes to be very inef- 
ficient. This is remarkably different from the unimolec- 
ular dissociation of more conventional, i.e. chemically- 
bounded, molecules, where statistical arguments are usu- 
ally applicable at least for the DOF directly involved 
in the fragmentation. In the classical description of 
the VP of Ne2Br2, the Br2 vibrational energy gradu- 
ally 'flows' into the vdW modes and eventually becomes 
large enough for dissociation to occur. A very limited 
number of the accessible final vibrational states is hence 
populated, the distributions being highly non-statistical 
(Fig. [3]). In addition, due to the very low density of 
vibrational states, the classical picture dramatically dif- 
fers from the quantum-mechanical description and hence 
the experimental observations. This is particularly clear 
from SB product branching ratios. In this case, the loss 
of 1.5 'quanta' provides enough energy to eject both Ne 
atoms up to v' = 22 (Table IIIII) , which is reflected in 
nonphysical trajectories dissociating via the Av' = — 1 
channel (Table ITV)) . The main advantage of GW over SB 
results in this case is to populate the qualitatively correct 
vibrational channels. However, as discussed above, the 
classical description is intrinsically inadequate and GW 
predicts a larger than observed propensity for Av' = — 1 
and —2 channels for the respective dissociation of the first 
and second Ne atoms. Also, the increasing importance of 
highly state-specific IVR (sparse regime) predicted in the 
experiment for states above v' = 19 cannot be adequately 
reproduced within QCT. Some additional discrepancies, 
e.g. in channels closings, are due to the PES employed 
and have been already discussed at the end of Sec. IIII CI 



2. Rotational distributions 

Together with vibrational distributions, product rota- 
tional state distributions provide a more detailed pic- 
ture of the VP process than dissociation lifetimes. In 
particular, they contain important information on the 
anisotropy of the PES, the vibrational level spacings and 
the effects of IVR and resonances on the fragmentation 
dynamics^*". Rotational distributions corresponding to 
various initial v' levels for the NeBr2 intermediate tri- 



13 



TABLE IV. Experimental'*^ and theoretical (this work) branching ratios in the vibrational predissociation of Ne2Br2(-B, u'). 



NeBr2(«' - n):NeBr2(i'' 


-1) 


Br2(«' — n): 


"QMv -2) 




v' W-2) (w'-3) {v' 


-4) 


{v'-l) W-i) 


(«'-4) C 


v' - 5) 


16 Experiment 




0.16 






QCT: SB/GW 0.04/0.27 




0.98/0.0 






17 Experiment 




0.31 


0.05 




QCT: SB/GW 0.08/0.29 




0.48/0.0 






18 Experiment 




0.33 


0.08 




QCT: SB/GW 0.16/0.35 




0.20/0.0 0.01/0.0 






19 Experiment 




0.45 


0.09 




QCT: SB/GW 0.28/0.45 




0.07/0.0 0.02/0.01 






20 Experiment 




0.65 


0.13 


0.03 


QCT: SB/GW 0.47/0.56 0.01/0.01 




0.02/0.0 0.05/0.03 






21 Experiment 




0.97 


0.29 




QCT: SB/GW 0.71/0.73 0.02/0.02 




0.01/0.0 0.10/0.06 






22 Experiment 




1.58 


0.51 


0.17 


QCT: SB/GW 1.04/0.92 0.06/0.04 




0.01/0.0 0.18/0.10 


0.01/0.0 




23 Experiment 




3.28 


1.37 


0.49 


QCT: SB/GW 1.52/1.22 0.15/0.09 0.01/0.0 


0.29/0.15 0.02/0.01 





atomic complex and the Br2 diatomic product are de- 
picted in panels a and b of Fig. [5] respectively. Only the 
GW results are shown because even if just a few vibra- 
tional states are populated, the SB and GW curves are 
only slightly different. This is not particularly surpris- 
ing; although the available energy more than halves in 
the range of v' considered (Table UlI)) . the tail of QCT 
distributions rapidly tend to zero before this effect be- 
comes important. 

The rotational state of the NeBr2 intermediate com- 
plex was calculated from j + ln2) = ~^2(i)5 which is valid 
only for J = 0. As seen in panel (a), the dissociation of 
the first Ne atom leaves the resultant triatomic complex 
in a highly-excited rotational state. The corresponding 
distribution extends over more than 30 rotational lev- 
els and peaks around 8-10. Such rotational excitation, 
together with the additional excitation of the vdW vi- 
brations, are expected to be the cause for the lifetime 
of the intermediate molecule being larger than that of 
the directly excited triatomic complex. This effect has 
been observed in the experiment^^ and is reproduced by 
our QCT calculations, as seen by comparing the corre- 
sponding lifetimes in Table ITTl with the QCT predictions 
in Table 1 of Ref. [H 

The GW rotational distributions for the Br2 product, 
shown in panel (b), extend over the full range of accessi- 
ble states determined by the maximum available energy 
in Table IIIIl For the sake of comparison, the initial dis- 
tribution P'^{j) has also been included in this panel. In 
the particular case of J = 0, symmetry considerations 
for the ground state are only responsible for even j val- 
ues contributing in (j)^. In general, the calculated 
rotational distributions are considerably hotter for the 



products, with a maximum at about j — 6-8 almost in- 
dependent of u', than for the initial states. The initial 
and final distributions have, nevertheless, similar shapes. 
Although experimental rotational distributions were not 
reported in the experiment^, a few points can be made 
with respect to QCT predictions: (a) the Boltzmann-like 
shape seems to be consistent with experimental measure- 
ments for the directly excited NeBr2 molecule^"*, and pre- 
vious QCT calculations^ii^; and (b) the Br2 product is 
more rotationally excited than in the photo-dissociation 
of NeBr^i^i^i^^, for which the distributions peak about 
f — 4-6. The extent of the QCT distributions may 
be expected to differ from future experimental measure- 
ments due to the PES employed (see discussion at the 
end of Sec. IIII Ct . particularly close to a channel closing. 

IV. SUMMARY AND CONCLUSIONS 

We have studied the vibrational predissociation (VP) 
of the Ne2Br2(-B, v' = 16-23) van der Waals (vdW) clus- 
ter using the quasi-classical trajectory (QCT) and the 
(Cartesian) coupled coherent states (CCS) methods. 

A sequential mechanism was used to fit the dynami- 
cal evolution of the different complexes involved. Both 
QCT and CCS are shown to provide very good esti- 
mates for the different dissociation lifetimes reported in 
the experiment"*^ and previous molecular dynamics with 
quantum transitions (MDQT) simulations^i. QCT pre- 
dictions are, however, obtained at a much lower com- 
putational expense. Various sources of non-exponential 
behavior have been identified and their implications ex- 
tensively discussed. In particular, the initial shape of 
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(a) NeBr2 intermediate 
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FIG. 5. (Color online) Rotational distributions in the vibra- 
tional predissociation of Ne2Br2(-B), for: (a) NeBr2 interme- 
diate complex; and (b) Br2 diatomic product. -P^(i) is the 
quantum distribution used to generate the quasi-classical ini- 
tial conditions. 



the QCT curves at short times arises from the 'non- 
democratic' selection of classical initial conditions and 
the classical description of the process. Meanwhile, 
the quantum CCS simulations give curves whose initial 
shape is in much better agreement with the experimental 
observations^''. However, the basis functions used to de- 
scribe the wave function follow dissociative trajectories. 
The high dimensionality of the phase space therefore re- 
sults in the basis functions decoupling and thus a semi- 
classical description of the cluster. This change is as least 
partially responsible for the departure at longer times of 
the curves from their initial exponential shape. An in- 
trinsic multi-step dissociation mechanism, as observed in 
the experiment-- and predicted by MDQT calculations-'^^, 
may also be responsible for such non-exponential behav- 
ior at long times. Additional calculations are needed in 
order to clarify this issue. The time scale over which the 
CCS gives a good description was increased significantly 
to around 15 ps by re-expanding the dissociated basis 
functions. Despite this improvement, the QCT simula- 
tions give curves that agree much more closely with the 
experimental results at longer times. The behavior of 
the QCT curves at longer times is, however, affected by 



low-dimensional tori in classical phase space. These tori 
are mainly in the form of quasi-periodic trajectories and 
are partially due to the weakness of vdW interactions. In 
both QCT and CCS calculations, the percentage of tra- 
jectories which correspond to the concerted mechanism 
is practically negligible, below 0.1% for QCT and about 
2% for CCCS. In this regard, it is important to note that 
trajectory-based approaches like CCS may become more 
effective than 'standard' quantum methods, mainly be- 
cause they allow analyzing the different mechanisms by 
simply inspecting the trajectories. 

As in our previous work with QCT on triatomic vdW 
molecules"^^, we found that application of the Gaussian 
weighting (GW) procedure yields survival probability 
curves, and consequently lifetimes, which are not signifi- 
cantly different from those calculated using the standard 
binning (SB) procedure. In addition, the capabilities of 
QCT in the description of the fragmentation kinetics was 
analyzed in detail by using reduced-dimensionality mod- 
els of the complexes and concepts from phase-space trans- 
port theory. 

We have reported QCT ro-vibrational product state 
distributions for the intermediate and final states of the 
VP process, and compared the vibrational distributions 
with the experimental results of Pio et al.'^'^ . As in pre- 
vious studies, e.g. Ref. [s^, the SB was found to popu- 
late nonphysical dissociation channels {Av' = — 1 in this 
case). This is due to energetic issues and can be easily 
solved by using the GW method. The latter, however, 
predicts a larger than observed propensity for dissoci- 
ation of the first (second) Ne atom via the Av' = — 1 
(-2) channel. We argued that this is a general problem 
in the classical description of the VP of vdW clusters, 
which may be attributed to the weakness of vdW inter- 
actions. Due to the latter, the characteristic frequen- 
cies for the diatomic subunit are usually one or more 
orders of magnitude larger than those of the bending 
mode within the vdW well. The energy transfer lead- 
ing to dissociation is inefficient and slow. To complicate 
matters further, as energy flows from the vibrational to 
the transitional modes, the vdW bending modes evolve 
from bounded, through hindered, to free type of motion 
and the role of many non-linear resonances become in- 
creasingly important^. This fact additionally worsens 
the quality of any classical description of the VP process, 
since quantum state-selectivity cannot be adequately de- 
scribed. As expected, the quantum-classical discrepan- 
cies will be more pronounced for relatively low v' , given 
the very low density of vibrational states. 

Dissociation of the first Ne atom leaves the NeBr2 in 
a highly-excited rotational state. This rotational excita- 
tion, and the additional excitation of the vdW vibration, 
are responsible for an increase dissociation lifetime of the 
intermediate complex, as compared to the lifetime of the 
directly excited NeBr2. This effect was observed in the 
experiment^ and is correctly predicted by our calcula- 
tions. Product Br2 rotational distributions are found to 
be significantly hotter in the products, extending over 
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the full range of available energies. 
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